arXiv: 1502.0621 lv4 [cond-mat.stat-mech] 31 Mar 2015 


Mortality, Redundancy, and Diversity in Stochastic Search 


Baruch Meerson^ and S. Redner^ 

^ Racah Institute of Physics, Hebrew University of Jerusalem, Jerusalem 91904, Israel 
^Department of Physics, Boston University, Boston, MA 02215, 

USA and Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico 87501, USA 

We investigate a stochastic search process in one dimension under the competing roles of mortality, 
redundancy, and diversity of the searchers. This picture represents a toy model for the fertilization 
of an oocyte by sperm. A population of N independent and mortal diffusing searchers all start at 
X = L and attempt to reach the target at a; = 0. When mortality is irrelevant, the search time 
scales as ro/lnA for InA ^ 1, where td ^ L^/D is the diffusive time scale. Conversely, when the 
mortality rate p of the searchers is sufficiently large, the search time scales as independent 

of N. When searchers have distinct and high mortalities, a subpopulation with a non-trivial optimal 
diffusivity are most likely to reach the target. We also discuss the effect of chemotaxis on the search 
time and its fluctuations. 
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Stochastic searching [l| underlies a wide variety of pro¬ 
cesses in bi olog y M , animal foraging , chemical re¬ 
actions |10l . [llll. and search operations for missing persons 
or lost items Umi. A basic goal is to minimize the time 
needed to successfully find a desired target. Through evo¬ 
lution and algorithmic developments, Nature and man 
have developed clever strategies to optimize searches. 
Many of these processes involve a single searcher that 
continues the search for as long as necessary to reach the 
target. In this Letter, we investigate the role of funda¬ 
mental influences on stochastic search that have not yet 
received sufficient attention: mortality, redundancy, and 
diversity. The notion that a searcher dies if it does not 
reach a target within a fixed time epitomizes unsuccess¬ 
ful foraging. While the properties of random walks that 
die at a fixed rate has been recently investigated [Tsl - fisll , 
the general problem of how to best conduct a search with 
mortal searchers is barely explored M- 

When a searcher can die before reaching a target, a 
natural way to promote success is by launching many 
searchers. This redundancy embodies the fertilization 
of an oocyte, where (in humans) roughly 3 x 10® sperm 
cells initially attempt to reach the oocyte after copula¬ 
tion. Fertilization is an extremely complex and multi¬ 
faceted process (see Refs. Sim for reviews). We make 
no pretense of accounting for the many steps that a sperm 
undergoes to reach and ultimately fertilize the oocyte. 
However, sperm mortality, redundancy, and diversity all 
play roles in this fertilization process. Inspired by these 
basic facts, we investigate an idealized scenario of fertil¬ 
ization that are driven by the above attributes. 

The geometry of the system is quite simple. An oocyte 
located at the origin, a; = 0, is represented by an ab¬ 
sorbing boundary condition. At time t = 0, N diffus¬ 
ing searchers, each of which represents a sperm cell, are 
launched from x = L. The basic questions that we ad¬ 
dress are: (i) What is the probability of fertilization— 
that at least one of the searchers reaches the origin? (ii) 
What is the average time to reach the origin as a function 
of N and the mortality rate of the sperm? If immortal 


searchers are uniformly distributed on the line, the sur¬ 
vival probability of the target in one dimension is known 
to decay as exp(—const, x y/t) (22l - E^ . However, the bi¬ 
ologically relevant situation where the searchers all start 
at the same point has been less extensively investigated 
even when the searchers are immortal Hi-Isl. 

It is worth noting two important points: First, if the 
searchers are immortal, then a basic result of first-passage 
processes is that a single searcher will eventually reach 
the target, but the average time for this successful search 
is infinite [s^. The average search time is aga.in infinite 
for = 2 searchers, but is finite for > 3 [^. Second, 
for mortal searchers that die at a fixed rate, the average 
time for a successful search, even for a single searcher, is 
finite and approaches zero as the mortality rate increases. 
In this high-mortality limit, the only way for a searcher 
to reach the target is to do so very quickly. 

We first determine how a single mortal diffusing 
searcher finds the target. Let p{x,t) be the probabil¬ 
ity density that the searcher is located at 0 < a; < oo at 
time t. Its dynamics is determined by the diffusion-decay 
equation dtp = Dd^p — p,p, where we assume a fixed dif¬ 
fusivity D and mortality rate p. For a target at a; = 0 
and a searcher starting at x = L, the image method 
gives the solution to this diffusion-decay equation as: 


p{x,t) = 


o-nt 


y/AirDt - 


^-(x-Lf/ADt _ ^-{x+Lf/ADt 
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The flux at a; = 0 yields the probability density that the 
target is reached at time t: 


fi{t) = D dxp{x,f)\x=o = 


L 


/ADt 

^ 7 
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while the probability Fi{t) = J* dt'fi{t') that the 
searcher reaches the target by time t is 


Fi{t) = i eVD ^ QYic{L j \/ ADt + y/pit) 

+ i e~'^^ erfc(L/V4£)t — y/Jjit) , (3) 
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where erfc(z) = 1 — erf(z), and erf z = e~'^^du 

is the error function. When the searcher is immortal 
(the case /r = 0), Eq. ([3|) reduces to the more familiar 
expression Fi (t) = erfc(L/\/4Dt) for the probability that 
a diffusing particle reaches the origin by time t. The 
probability = Fi{t = oo) that the target is eventually 
found is 





dtfi{t) 


= e 


y/JI/DL 


( 4 ) 


For an immortal searcher = 1, whereas for high mor¬ 
tality rate Fi is exponentially small. 

From the first-passage distribution ([2]) , the average and 
the variance of the search time are 


(T) 


J^dtfiit) 
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where to = L'^ j AD is the diffusion time and = l//i 
is the lifetime of a mortal searcher. Notice that low- 
order moments of the search-time distribution involve 
non-trivial combinations of to and r^. Also notice that 
as /r —>■ 0, corresponding to immortal searchers, both (T) 
and E^ diverge. 

We now turn to the case of N immortal searchers. The 
probability that one of them first reaches the target at 
time t is 


fN{t) = Nh{t)[l-F^{t)] 
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where the third line follows from the z -C 1 asymptotic 
of erf(z). An important feature of /Ar(t) is that it has 
the algebraic long-time tail /jv(t) ^ Thus, as 

mentioned in the introduction, the average search time 
is divergent for < 2, but finite for > 3 [^ . 

We now determine the average time for the first out 
of A^ > 3 immortal searchers to reach the target. By 
definition 
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We i nteg rate by parts and define the scaled variable z = 
Ll\/ADt to recast this expression as [s^ 

/■oo 

{Tn)=toA>{N), A>{N) = 2 — erf^z. (8) 

Jo ^ 


For very large N, erf^z effectively becomes the Heavi¬ 
side step function 0(z —Zq), where Zq = zo{N) :$> 1. To 
determine zo{N), we use the large-z asymptotic erf z ~ 
1 — e“^^/(i/^z), so that erf’^z ~ exp[—A^e“^^/(v^-z)]- 

Now zo{N) is determined, with logarithmic accuracy, 

2 2 

from the condition Ne~^ /z ^ 1, or ze^ ^ A^. This 
yields zq ~ [VF(2A^^)]^/^/-v^, where W{-) is the Lambert 
W-function—the inverse of f{W) = We^ [s^. Thus 


A>{N )^2 r ^ 
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To leading order in InA^ ^ 1, this yields 


(Tn) 


AD In N ■ 


( 10 ) 


This result coincides with the asymptotic that was 
quoted in Ref. [s^ • The average search time (flUl) decays 
only logarithmically with N but still gives to a reduction 
by a factor of 20 compared with the characteristic diffu¬ 
sion time To = /AD for the typical number of human 
sperm {N = 3 x 10®) that attempt to fertilize an oocyte. 

For N identical and mortal searchers that all start from 
X = L, the probability that at least one of them eventu¬ 
ally reaches the target is 


p = Yl = 1 - (1 - ( 11 ) 

k=l ^ 7 


with Fi given by Eq. ((T|) . It is convenient to introduce the 
dimensionless mortality rate M = ^JJPro. For high mor¬ 
tality, M 1, we approximate p ~ 1 — exp(— 
which changes rapidly from being vanishingly small for 
N < Nc{M) to being close to 1 for A > Nc{M), with 
Nc{M) ~ ^ 1. An important message from this 

simple argument is that a huge redundancy of searchers 
is needed to offset their high mortality for a search to be 
successful. 

For mortal searchers, reaching the target is not guar¬ 
anteed, and the average time of successful search is 


. _ Io^tMt)dt ^ vhi(A,Af) 


( 12 ) 


Here «'fe(A,M) = dz z-^^ and, by using 

Eqs. dS]) and (O, z = L/y/ADt, and also expressing all 
variables in scaled form, straightforward algebra gives 


^{N,M,z) = z^-h 
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We are interested in the high-mortality regime, M ^ 1, 
where the probability that a single searcher eventually 
reaches the target, as given in Eq. ®, is exponentially 
small. The inequality M ^ 1 enables us to evaluate the 
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integral in by the standard Laplace method [s^. The 
saddle point z^{N, M) is found by minimizing $(iV, M, z) 
with respect to z: 

dm,M,z) 2PP 2(iV-l) 

dz “ z3 ® ^ ^ 


where we again use the z ^ 1 asymptotic of erfc( 2 :). 
Since the functions 4'o and 'I'l include the same exponent 
g$(Af,M.z)^ we obtain, after cancellations, d'l/'I'o — z~^. 

There are two distinct limiting behaviors for the saddle 
point z* that depend on the interplay between N and 
M. For N <g; , z^^ is determined by balancing 

the hrst two terms in Eq. (ITSl) . This yields z* ~ y/M. In 
this regime, we obtain 



(a) (b) 


FIG. 2: (Color online) (a) Non-monotonic dependence of {Tn) 
from Eq. in on M for N = 10. The large-M asymptotic 
{TN)/rD zs 1/M (Eq. (fid)) ') is shown dashed, (b) {Tn)/td 
vs. InA^ for M — 2. The dashed and dotted lines are the 
respective asymptotics (0) and CQl). 


{Tn) ^ y/Tf,TD 


Tp 

M 


L 
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which coincides with Eq. (15a)) for the average search time 
of a single searcher. In this high mortality regime, only 
the fastest searcher contributes to the average time, while 
the rest of the searchers (unless their number is huge, see 
below) are superfluous. 

In the limit of y/Me'^^, but with M also large, the 
saddle point is given by 2z — 2{N— l)e~^ /ypK ~ 0, which 
yields z* ~ zg, the same quantity that arises for immor¬ 
tal searchers, as given above Eq. As a result, the 
final expression for (Tjv) coincides with Eq. (flOl) . Thus 
when the number of searchers is extremely large, their 
mortality (even when relatively high) is irrelevant in the 
determination of average search time. 
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FIG. 1: (Color online) Contour plot of {Tn)/td from Eq. (1121) 
vs. M and In A. 

Figure [D shows the dependence of (Tn) on M and In N. 
For not too large N, (Tn) in Eq. (1121) exhibits a maxi¬ 
mum as a function of the mortality rate, as evident from 
Fig. [T] and Fig. [IJa). This latter figure also shows that 
the M 1 asymptotic {Tn) — rp/M is accurate already 
at M ~ 2. Figure [2Kb) shows {Tn)/tp vs . In A and the 
asymptotics (O and (flUl) . 


The biologically relevant situation where mortal 
searchers have distinct diffusivities Dk, k = 1,..., N [s^ 
gives rise to a new optimization criterion for a successful 
search. Now the probability that a searcher with diffu- 
sivity Dk eventually reaches the target is 

D{Dk) = < 1. (15) 


When NiF{Dk) <C 1 for typical diffusivity values, we 
may approximate the probability that a single searcher 
reaches the target as p ~ T{Dk). For A 1, we 

may make the simplifying assumption that the diffusiv¬ 
ity distribution is a Gaussian that is centered about a 
typical value Dq. P{D) = ^ ^-{D-Pa) jia ^ 

a Dq. Then the summation for p may be replaced by 
the integral 


p~ / dDT{D)P{D) 

Jo 
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We again evaluate the integral by the Laplace 
method [36l| . By straightforward rescaling, we redefine 
the above integral as 



(17) 


where $(z) = z 3/2 _|_ with 


A = 


^0 

2a2 


y/Pp 

Ly/Ji 


and z 
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We find the saddle point z* from the equation d>'(z) = 
2A(z — 1) — (1/2 )z“ 3/2 = 0. The exact solution of this 
equation is cumbersome, and we confine ourselves to the 
asymptotics in the limits A ;g> 1 and A <C 1. In the low 
mortality limit where A 1 (but still L^J p/Ao 1), z* 
is very close to 1, and we find z* ~ 1 -|- (4A)“3. Applying 
the Laplace method, the arrival probability is 


p ~ exp 



A» 1. 


(18) 
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In this case, the dominant contribution to p comes from 
“typical” searchers—those whose diffusivity is close to 
Dq. This behavior is to be expected, as in the low- 
mortality limit the Gaussian in the integrand of Eq. (USD 
effectively acts as a delta-function peak centered at Dq. 

The high-mortality limit, A <C 1, is more interesting. 
Here z* ^ 1, and we can replace z — 1 with z in 4>(z) to 
arrive at z* ~ In this case, a small number of 

highly active searchers with D ^ Dq give the dominant 
contribution to the probability that the target is found. 
Performing the Gaussian integral, we obtain 


p ~ 



(21/5 + 2-9/5) 


A<1. 


... 

Surprisingly, this result for the arrival probability is in¬ 
dependent of the average diffusivity Dq. Numerical inte¬ 
gration of m shows that the arrival probability decays 
monotonically with the basic dimensionless parameter A, 
and the asymptotic forms (IT^ and (US match the full 
solution in the respective limits of H I and H ^ I. 

There are two additional attributes that naturally arise 
in the phenomenon of oocyte fertilization by sperm that 
can be accounted for within the present framework. One 
is the possibility of a diversity in searcher lifetimes. If 
the mortality rate p is normally distributed, we find that 
the search dynamics is only trivially affected [s^. In 
this case, the main contribution to the probability that 
an ensemble of searchers reaches the target comes from 
typical searchers. 

The role of chemotaxis, which is known to be a dom¬ 
inant effect as sperm approach the fertilization site [ 2 ^, 
has more interesting consequences. A simplistic way to 
model chemotaxis is to include a constant drift velocity v 
in the diffusion-decay model. The probability evolution 
of a single searcher is now governed by the equation of 
motion dtp — vdxP = —pp + Dd^p, where u > 0. When 
a single searcher starts at x = L and seeks a target that 
is at the origin, the solution to this equation, subject to 
the boundary conditions p{x = 0,t) = Pi^ = 00 , t) = 0 
and the initial condition p{x, t = 0) = S{x — L), can again 
be found by the image method |34| . The result is 


p{x,t) = 


\/A'KDt 




.^ — {x+vt—L)‘^/ADt 


_^vL/D ^-(x-i-vt+L)^/4Dt 


( 20 ) 


From this expression, the probability density that the 
target is reached at time t is 


hit) 


'J AttDD 


( 21 ) 


Gonsequently, the probability that the target is eventu¬ 
ally reached is 


T, = r dt hit) = e-L/2D ^22) 
Jo 

From Eq. (ED, the average search time and its variance 
are 


(T) = , , S" =-—. (23) 

Interestingly, the variance is non-monotonic in I?, and 
the largest fluctuations in the search time occur when 
D = z;^/(2/i). Using our general approach, we can ex¬ 
tend this chemotactic search process to the case of many 
searchers and determine the behavior of the search time 
on the N, M, and the rescaled drift velocity, the Peclet 
number [^. It is also possible to investigate this search 
process in higher dimensions and in more realistic ge¬ 
ometries, with the goal of providing a realistic but still 
tractable model for fertilization. 

To summarize, we elucidated the competing roles of 
mortality, redundancy, and diversity on a search pro¬ 
cess that represents a caricature for the fertilization of 
an oocyte by sperm. To optimize this search, one strat¬ 
egy is to have a sophisticated search algorithm. However, 
nature often seems to prefer the brute-force approach 
of dispatching many almost-identical searchers that fol¬ 
low a simple search algorithm. The effectiveness of this 
redundancy is counterbalanced by the mortality of the 
searchers, and it is miraculous that the correct (and a 
very small) number of sperm actually reach the oocyte 
in human fertilization. While we do not offer insight 
into why this miracle occurs, we quantified the dynam¬ 
ics of this search process as a function of the number of 
searchers and their mortality rate. We also found that 
searcher diversity can compete with mortality so that 
only the most active searchers are successful. 
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